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Impacts  on  planetary  bodies  can  lead  to  both  prompt  secondary  craters  and  projectiles  that  reimpact  the 
target  body  or  nearby  companions  after  an  extended  period,  producing  so-called  “sesquinary”  craters. 
Here  we  examine  sesquinary  cratering  on  the  moons  of  Mars.  We  model  the  impact  that  formed 
Voltaire,  the  largest  crater  on  the  surface  of  Deimos,  and  explore  the  orbital  evolution  of  resulting 
high-velocity  ejecta  across  500  years  using  four-body  physics  and  particle  tracking. 

The  bulk  of  mass  transfer  to  Phobos  occurs  in  the  first  102  years  after  impact,  while  reaccretion  of  ejecta 
to  Deimos  is  predicted  to  continue  out  to  a  104year  timescale  (cf.  Soter,  S.  [1971].  Studies  of  the 
Terrestrial  Planets.  Cornell  University).  Relative  orbital  geometry  between  Phobos  and  Deimos  plays  a 
significant  role;  depending  on  the  relative  true  longitude,  mass  transfer  between  the  moons  can  change 
by  a  factor  of  five.  Of  the  ejecta  with  a  velocity  range  capable  of  reaching  Phobos,  25-42%  by  mass  reac¬ 
cretes  to  Deimos  and  12-21%  impacts  Phobos.  Ejecta  mass  transferred  to  Mars  is  <10%. 

We  find  that  the  characteristic  impact  velocity  of  sesquinaries  on  Deimos  is  an  order  of  magnitude 
smaller  than  those  of  background  (heliocentric)  hypervelocity  impactors  and  will  likely  result  in  different 
crater  morphologies.  The  time-averaged  flux  of  Deimos  material  to  Phobos  can  be  as  high  as  11%  of  the 
background  (heliocentric)  direct-to-Phobos  impactor  flux.  This  relatively  minor  contribution  suggests 
that  spectrally  red  terrain  on  Phobos  (Murchie,  S.,  Erard,  S.  [1996].  Icarus  123,  63-86)  is  not  caused  by 
Deimos  material.  However  the  high-velocity  ejecta  mass  reaccreted  to  Deimos  from  a  Voltaire-sized 
impact  is  comparable  to  the  expected  background  mass  accumulated  on  Deimos  between  Voltaire-size 
events.  Considering  that  the  high-velocity  ejecta  contains  only  0.5%  of  the  total  mass  sent  into  orbit,  ses¬ 
quinary  ejecta  from  a  Voltaire-sized  impact  could  feasibly  resurface  large  parts  of  the  Moon,  erasing  the 
previous  geological  record.  Dating  the  surface  of  Deimos  may  be  more  challenging  than  previously 
suspected. 

Published  by  Elsevier  Inc. 


1.  Introduction 

Several  features  about  the  surface  geology  on  the  moons  of 
Mars  remain  poorly  understood.  The  grooves  on  Phobos,  which 
do  not  exist  on  Deimos,  have  received  the  most  attention 
(Horstman  and  Melosh,  1989;  Thomas,  1979;  Weidenschilling, 
1979),  and  theories  for  their  formation  continue  to  be  proposed 
to  this  day  (Asphaug  et  al.,  2015b;  Basilevsky  et  al.,  2014; 
Hamelin,  2011;  Murray  et  al,  2006;  Wilson  and  Head,  2015; 
Nayak  and  Asphaug,  2015).  However  this  is  far  from  the  only  mys¬ 
tery.  Though  both  moons  are  heavily  cratered,  with  saturated 
surfaces  and  fine-grained  regolith  from  impact  debris  accumula¬ 
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tion  (Lunine  et  al.,  1982;  Thomas,  1979),  a  large  portion  of  ejecta 
produced  on  Deimos  is  retained  in  the  form  of  crater  fill  of  ~5  m 
depth,  a  phenomenon  not  noted  on  Phobos  (Thomas  and 
Veverka,  1980a).  This  difference  is  still  unexplained  (Lee,  2009). 
The  surface  of  Deimos  is  also  significantly  smoother  and  brighter 
than  Phobos,  likely  a  result  of  crater  fill  (Thomas,  1993;  Thomas 
et  al.,  1996;  Veverka,  1978). 

Phobos  also  exhibits  two  distinct  spectral  units,  one  of  “redder” 
origin  and  one  of  “bluer”  origin,  possibly  stemming  from  a  compo¬ 
sitional  difference  (Lee,  2009;  Murchie  and  Erard,  1996).  The  bluer 
unit  is  associated  with  the  Stickney  crater  and  an  origin  from 
depth.  The  redder  unit  associated  with  the  surface  and  small  cra¬ 
ters;  it  is  spectrally  similar  to  D-type  asteroids,  but  also  to  Deimos 
(Murchie  and  Erard,  1996, 1993).  It  has  been  proposed  that  the  red 
unit  is  a  wide-spread  shallow  layer  superimposed  on  a  blue  base 
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(Murchie  and  Erard,  1996),  for  which  there  are  four  possible  causes 
(Britt  and  Pieters,  1988;  Murchie  et  al.,  1991):  (1)  accretion  of 
D-asteroid  material  onto  blue  Phobos  material;  (2)  optical  alter¬ 
ation  of  the  bluer  unit;  (3)  accretion  of  ejecta  from  Martian  basin 
impacts  and  subsequent  space  weathering  or  (4)  Phobos  is  an 
inherently  heterogeneous  rubble  pile  and  the  red/blue  units  are 
end-member  compositions.  One  aim  of  our  study  is  to  investigate 
the  possibility  that  the  red  veneer  on  top  of  the  base  blue  unit 
may  be  ejecta  accreted  from  Deimos  rather  than  Mars. 

Previous  work  has  established  that  impact  ejecta  can  reimpact 
the  target  body  or  nearby  companions  after  an  extended  period,  cre¬ 
ating  so-called  “sesquinary”  impact  morphology.  Examples  of  ses- 
quinary  studies  in  the  literature  include  Earth’s  Moon  (Gladman 
et  al.,  1995),  Phobos  (Nayak  and  Asphaug,  2015),  Io  (Alvarellos 
et  al.,  2008),  Ganymede  (Alvarellos  et  al.,  2002),  Europa  (Zahnle 
et  al.,  2008)  and  Pluto  (Bierhaus  and  Dones,  2014).  For  Mars,  previ¬ 
ous  work  suggests  ejecta  released  at  slightly  greater  than  the  satel¬ 
lite’s  escape  velocity  could  remain  in  the  system  and  subsequently 
reimpact  at  low  relative  velocities  (Soter,  1972, 1971).  Possible  evi¬ 
dence  for  this  was  noted  in  analysis  of  Viking  images  (Veverka  and 
Duxbury,  1977),  however  the  efficiency  of  this  process  was  previ¬ 
ously  unknown  (Thomas,  1979).  We  report  here  on  the  distribution 
of  impact  velocities  and  geometries  from  inter-moon  mass  transfer 
trajectories,  and  present  conclusions  on  the  role  and  importance  of 
sesquinary  mass  transfer  between  the  Martian  moons. 


a  = 


gok 


i2Z+l 


4Z(Z  -  2) 


(4) 


Using  a  final  crater  radius  R/=  1500  m  for  Voltaire  (Veverka, 
1978),  the  transient  crater  radius  is  calculated  as  Rt  =  0.65  Rf 
(Barnhart  and  Nimmo,  2011).  For  the  analysis  presented  here  the 
number  of  streamlines  (n)  has  been  chosen  to  yield  a  suitably 
dense  streamline  distribution  with  velocities  greater  than  the  Dei¬ 
mos  escape  velocity.  Setting  Rmin  =  0  and  varying  Rmin  ^  r  ^  Rt, 
n  =  600  streamlines  evenly  spaced  in  radius  are  generated  within 
the  Voltaire  crater.  Converting  streamlines  into  axisymmetric  coor¬ 
dinates  (cf.  Barnhart  and  Nimmo,  201 1 ,  Fig  1 ),  we  extract  the  radial 
and  vertical  coordinates  as: 


r  =  Ri  sin  0(1  -cos0)z-2 


(5) 


z  =  Ri  cos  0(1  -  cosO)z~2 

where  9  is  the  angle  from  the  vertical  {6\6  e  0  :  n/2}  and: 

n  ~  Rmin 


(6) 

(7) 


2.2.  Creating  3-D  velocity  streamlines 


2.  Methods 

2.2.  Impact  model:  generating  2-D  velocity  streamlines 

Voltaire,  the  largest  confirmed  crater  on  Deimos,  has  a  diameter 
of  3  km  (Thomas  and  Veverka,  1980b;  Thomas,  1979).  By  modeling 
the  orbital  evolution  of  ejecta  from  the  Voltaire-forming  impact, 
we  aim  to  characterize  an  end-member  case  of  mass  transfer  from 
Deimos  to  other  Martian  system  bodies. 

To  model  the  streamlines  ejected  by  the  Voltaire  impact,  we  use 
a  simplified  form  of  Maxwell’s  Z-model  (Maxwell  and  Seifert, 
1974;  Maxwell,  1977;  Roddy,  1977).  Though  limited  by  its  neglect 
of  interactions  across  streamlines,  the  Z-model  reasonably  approx¬ 
imates  several  experimentally  observed  features  (Melosh,  1989; 
Richardson  et  al.,  2007).  The  limitations  of  a  Z-model  implementa¬ 
tion  are  discussed  at  length  by  (Barnhart  and  Nimmo,  2011).  Our 
application  is  only  concerned  with  ejecta  streamlines  that  escape 
Deimos,  and  is  unaffected  by  the  details  of  cratering  flow  beneath 
the  ground  plane,  surface  material  mixing  during  ejection  or  direct 
retention  and  emplacement  of  deposits.  Therefore  it  provides  a 
suitable  level  of  insight  into  an  outbound  velocity  distribution; 
approximations  made  by  the  Z-model  are  unlikely  to  alter  our 
qualitative  results. 

We  adopt  the  formulation  of  Barnhart  and  Nimmo  (2011),  who 
use  Z  =  2.71  for  a  Mars  application.  When  tested  against  numerical 
computations,  Z  =  2.7  represents  surface  explosion  cratering  flow 
well  (Melosh,  1989).  All  streamlines  are  ejected  at  a  constant  angle 
of  35.4°  from  the  horizontal,  set  according  to  the  relation  (Maxwell, 
1977): 

sej  =  tan1  (Z  -  2)  (1) 

Outbound  radial  ( vr )  and  vertical  ( vz )  ejection  velocities  vary 
inversely  with  distance  from  the  center  of  the  crater  r  (Maxwell, 
1977)  as: 

vr  =  a/rz  (2) 

vz  =  (Z-2)vr  (3) 

where  gD  denotes  the  acceleration  due  to  gravity  for  Deimos 
(0.003  m/s2)  and: 


The  2-D  axisymmetrical  distribution  is  now  used  to  create  an 
approximation  to  a  3-D  excavation.  The  fate  of  the  ejecta  particle 
(reaccretion  to  Deimos,  impact  to  Phobos,  impact  to  Mars  or 
escape)  can  vary  greatly  depending  on  the  azimuth  of  the  stream¬ 
line.  To  rotate  around  the  azimuthal  direction,  we  define  the 
Topocentric  Horizon  frame  (Appendix  A),  adapted  from  the 
South-East-Zenith  (SEZ)  frame  (Vallado,  2013).  The  azimuth  of 
the  ejection  velocity  vector  p  is  measured  from  the  north,  clock¬ 
wise  as  viewed  from  above  the  impact  site.  We  select  a  30°  span 
as  a  compromise  between  computational  efficiency  and  sampling 
a  variety  of  azimuths  across  the  possible  solution  space,  such  that 
P\P  e  (0  :  7i/6  :  2n)  for  a  total  of  11  possible  azimuths.  This  yields  a 
three-dimensional  outbound  velocity  distribution  tied  to  Voltaire. 
For  use  with  the  Mars  gravity  system  integrator,  these  coordinates 
are  then  rotated  into  the  Mars  Centered  Inertial  (MCI)  frame; 
details  of  coordinate  transformations  through  the  Deimos- 
Centered  Deimos-Fixed  (DCDF)  and  Deimos-Centered  Inertial 
(DCI)  frames  are  presented  in  Appendix  A. 

Finally,  we  are  specifically  interested  in  those  streamlines  that 
have  sufficient  velocity  to  reach  the  orbit  of  Phobos.  Since  both 
moons  lie  in  the  same  orbital  plane  (Cazenave  et  al.,  1980),  the 
minimum  velocity  at  Deimos  to  reach  Phobos  can  be  analytically 
calculated  with  the  Hohmann  transfer  (Section  6.3,  Curtis,  2013). 
Particles  begin  to  cross  the  orbit  of  Phobos  at  velocities  above 
500  m/s,  so  we  set  the  lower  bound  on  velocities  of  interest  at 
400  m/s.  From  Deimos,  the  minimum  velocity  to  escape  the  gravi¬ 
tational  well  of  Mars  is  analytically  approximated  as  (Eq.  (2.80), 
Curtis,  2013): 

Vesc  —  Deimos  (8) 

where  rDeimos  is  the  distance  from  Deimos  to  Mars  and  p  is  the  pro¬ 
duct  of  the  gravitational  constant  and  the  mass  of  Mars.  From  (8), 
vesc  =  1.91  km/s;  we  set  the  upper  bound  on  velocities  of  interest 
at  2  km/s.  Therefore,  we  examine  velocity  streamlines  in  the 
range  {v\ v  e  400  :  2000  m/s}.  Nineteen  of  600  streamlines  fall 
within  this  range;  rotated  around  11  azimuthal  positions,  this 
creates  a  209-streamline  distribution.  While  we  focus  here  on 
ejecta  with  sufficient  velocity  to  reach  Phobos  (~400  m/s),  we  note 
that  the  majority  of  ejecta  launched  from  Deimos  at  lesser  velocities 
will  ultimately  re-impact  Deimos. 
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2.3.  Planetary  system  model 

Next  we  detail  the  formulation  of  the  planetary  model.  Cen¬ 
tered  at  the  primary,  the  Mars  gravity  system  is  modeled  with 
12x12  gravity  harmonics  from  the  NASA  Planetary  Data  System 
[pds-geosciences.wustl.edu]  (Murchie,  2010).  The  effects  of  perma¬ 
nent  solid  tides  are  included,  truncated  to  the  size  of  the  gravity 
field.  The  present-day  orbit  of  Deimos  is  likely  similar  to  its  pri¬ 
mordial  orbit  (Burns,  1978;  Lambeck,  1979);  we  generate  the 
500-year  orbits  of  Phobos  and  Deimos  analytically  from  modern- 
day  mean  orbital  parameters  [ssd.jpl.nasa.gov,  Table  1].  A  subset 
of  streamlines  was  run  against  high-precision  orbits  generated 
for  Phobos  and  Deimos  [Genova  and  Folkner,  personal  communica¬ 
tion,  2015];  results  were  not  found  to  differ  substantially  from 
those  run  against  the  analytical  orbits.  In  the  interest  of  computa¬ 
tional  speed,  we  adopt  the  analytical  formulation  hereafter. 

Due  to  its  proximity  to  Mars,  the  orbit  of  Deimos  is  primarily 
influenced  by  Mars’  oblateness;  the  third-body  effect  from  the 
Sun  or  other  planetary  bodies  such  as  Jupiter  is  negligible  (Burns, 
1972).  Ejecta  released  from  the  orbits  of  Deimos  will  follow  a  sim¬ 
ilar  pattern;  we  therefore  neglect  these  third-body  effects.  Simi¬ 
larly,  solar  radiation  pressure  is  a  second-order  effect  when 
compared  to  solar  gravity  perturbations  (Farnocchia  et  al.,  2014; 
Klacka,  2002);  we  neglect  this  effect  as  well.  However  for  complete 
understanding  of  orbital  evolution  within  the  Martian  system  we 
include  third-body  perturbation  effects  from  Phobos  and  Deimos, 
making  the  physics  of  our  model  a  four-body  problem. 

We  assume  that  any  particle  that  enters  the  Hill  sphere  (Table  1 ) 
of  either  Moon  will  be  captured  by  it.  Due  to  its  irregularly  triaxial 
shape,  Deimos  has  an  uneven  gravity  field  that  causes  the  escape 
velocity  to  be  lower  at  the  sub-Mars  and  anti-Mars  points  (Davis 
et  al.,  1981).  Ejecta  in  the  4-6  m/s  range  will  see  the  largest  variation 
in  range  (Thomas,  1993);  since  the  slowest  particle  we  consider  is 
ejected  at  ^400  m/s  and  the  escape  velocity  varies  on  the  order  of 
cm/s,  we  can  safely  assume  that  the  escape  velocity  at  Voltaire 
equals  the  average  escape  speed  over  Deimos. 

Finally  we  consider  the  impact  of  relative  orbital  geometry. 
Though  an  analytical  formulation  has  been  used  to  consider  similar 
problems  in  the  past  (Dobrovolskis  and  Burns,  1980;  Soter,  1971; 
Thomas,  1998),  this  approach  may  be  insufficient  for  a  full  under¬ 
standing  of  ejecta  dynamics.  Phobos  is  closer  to  Mars  than  any 
other  planetary  satellite,  and  is  the  only  Moon  with  an  orbital  per¬ 
iod  less  than  the  rotational  period  of  its  primary  body  (Burns, 
1972).  The  flux  of  material  impacting  Phobos  can  vary  drastically 
between  inferior  and  superior  conjunctions  between  Phobos  and 
Deimos.  The  difference  in  the  true  longitude  between  Deimos 
and  Phobos  can  (and  does,  see  Fig.  1 )  change  the  outcome  of  a  Pho¬ 
bos  collision  to  a  Mars  collision,  or  vice  versa.  Therefore,  though 
the  orbits  of  Phobos  and  Deimos  are  generated  analytically,  all 
propagation  in  this  work  is  ephemeris-centered. 

Deimos  has  an  orbital  period  of  30.3  h,  and  Phobos  7.5  h.  To 
evaluate  the  fate  of  ejecta  across  the  range  of  possible  Mars- 
Deimos  orbital  geometries,  we  discretize  this  orbit  into  28  geome¬ 
try  configurations  (GCs),  evenly  spaced  in  one-hour  increments 
from  the  Deimos  apoapsis.  In  this  time  period,  Phobos  completes 
nearly  four  orbits  of  Mars,  allowing  for  discretization  of  the  range 
of  possible  Mars-Phobos-Deimos  orbital  geometries  as  well.  At 
each  GC,  209  streamlines  are  released  and  propagated,  for  a  total 
of  5852  streamlines,  thereby  ensuring  a  robust  capture  of  the 
impacting  process  despite  variations  in  orbital  positions  and 
conjunction  geometries. 

1  To  test  the  validity  of  using  the  Hill  sphere  as  an  impact  boundary,  we  selected 

100  Phobos  impact  trajectories  at  random  and  integrated  them  with  an  impact 

boundary  of  13  km,  the  longest  semi-major  axis  of  Phobos  (Murchie  et  al.,  2003); 
98/100  trajectories  were  still  found  to  impact.  Therefore,  our  qualitative  results  are 

not  impacted  by  the  use  of  the  Hill  sphere  as  an  impact  boundary. 


Table  1 

Mean  orbital  parameters  and  constants  for  Phobos  and  Deimos. 


Deimos 

Phobos 

Semi-major  axis 

23,485  km 

9389.8  km 

Eccentricity 

0.00115571 

0.0164255 

Inclination 

1.79  deg 

1.09  deg 

Right  ascension  of  ascending  node 

148.0  deg 

319.9  deg 

Argument  of  periapsis 

123.3  deg 

270.7  deg 

Mean  longitude3 

109.7  deg 

190.6  deg 

Rate  of  mean  longitude 

0.00330049  deg/s 

0.0130317  deg/s 

Mean  radius 

6.2  km 

1 1 .3  km 

Acceleration  due  to  gravity 

0.003  m/s1 2 

0.0057  m/s2 

Escape  velocity 

5.56  m/s 

11.39  m/s 

Hill  sphere  radius 

16.5  km 

21.5  km 

a  Mean  longitude  is  calculated  with  reference  to  a  mean  epoch  coordinate  sys¬ 
tem:  The  mean  equator-mean  equinox  coordinate  system  is  evaluated  at  the  epoch 
of  the  object.  The  starting  epoch  is  arbitrary  due  to  our  evaluation  of  the  orbital 
dynamics  at  multiple  relative  geometry  configurations  between  Mars,  Phobos  and 
Deimos  that  encompass  all  possible  geometries  between  the  bodies. 

Estimates  for  lifetime  of  ejecta  in  the  Martian  system  range 
from  102  to  104  years  (Davis  et  al.,  1981;  Soter,  1971).  We 

propagate  each  streamline  in  the  Mars  gravity  system  for  tmax  = 
500  years,  stopping  sooner  only  in  the  event  of  a  planetary  body 
collision  or  departure  from  the  Mars  gravitational  sphere  of  influ¬ 
ence.  This  length  of  integration  balances  computational  feasibility 
with  permitting  a  statistically  significant  number  of  ejected  parti¬ 
cles  to  impact  or  escape.  As  we  shall  show,  the  uncertainty  intro¬ 
duced  by  doing  so  does  not  affect  our  conclusions.  It  also  permits 
the  use  of  a  Runge-Kutta  integrator  without  excessive  approxima¬ 
tions  to  the  perturbed  Hamiltonian  (Leimkuhler  and  Reich,  2004). 
A  seventh-order  Runge-Kutta  integrator  with  eighth-order  error 
control  is  used  for  all  orbit  propagations.  The  maximum  permitted 
relative  error  is  10-10.  The  Tisserand  parameter  is  used  to  evaluate 
the  performance  of  the  integrator,  according  to  which: 

\+2^  -  e2)  cos  i  2*  constant  (9) 

where  a,  e  and  i  are  the  semi-major  axis,  eccentricity  and  inclination 
of  the  orbit.  The  differences  in  the  Tisserand  parameter  for  an  indi¬ 
vidual  particle  are  no  greater  than  10-5,  i.e.,  at  most,  a  0.001% 
change  across  the  integration  period. 

3.  Results 

Fig.  1  plots  the  fate  of  5852  massless  ejecta  particles  in  the  Mar¬ 
tian  gravity  system  across  500  years  as  a  function  of  the  orbital 
geometry  and  Voltaire  ejection  velocity.  The  number  of  reaccre¬ 
tions  to  Deimos  and  particles  still  flying  are  relatively  constant, 
with  minor  fluctuations.  However,  if  the  Voltaire  impact  occurs 
when  Deimos  is  near  periapsis,  a  spike  in  Mars  impacts  and  a  cor¬ 
responding  drop  in  particles  escaping  the  system  are  noted.  This  is, 
in  fact,  due  to  the  difference  in  true  longitude  between  Phobos  and 
Deimos  at  the  time  of  ejecta  launch;  depending  on  the  relative  con¬ 
junctive  geometry,  the  mass  flux  from  Deimos  to  Phobos  can  be  up 
to  500%  higher. 

It  is  surprising  to  note  that  there  is  more  mass  flux  from  Deimos 
to  Phobos  as  opposed  to  Mars;  intuitively,  one  would  expect  that 
most  mass  ejected  from  Deimos  would  either  reaccrete  or  spiral 
down  to  Mars.  Because  greater  mass  is  released  at  lower  ejection 
velocities  (O’Keefe  and  Ahrens,  1985),  these  results  suggest 
impacts  on  Deimos  may  have  an  effect  on  Phobos’  geology;  we 
shall  attempt  to  estimate  the  magnitude  of  that  effect  in  Section  4. 

For  impacts  with  Mars,  Phobos  or  Deimos,  impact  velocity  is 
calculated  with  reference  to  the  Planet-Centered  Planet  Fixed 
frame  in  question  (Appendix  A).  We  chart  the  variation  of  impact 
speeds  at  Phobos  and  Deimos  across  28  GCs  (Fig.  2).  For  Phobos, 
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Fate  of  ejecta:  Function  of  orbital  position 
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Fig.  1.  The  fate  of  Deimos  ejecta  from  Voltaire  as  a  function  of  (Left)  Orbital  Geometry  and  (Right)  Ejection  Velocity  from  Voltaire.  GC-1  is  near-apoapsis;  GC-14  is  near- 
periapsis.  Note  that  the  difference  in  the  Deimos-Phobos  angle  is  the  important  quantity  with  regard  to  ejecta  fate  (see  text). 


regardless  of  orbital  location  at  the  time  of  ejecta  release,  impact 
velocity  scales  linearly  with  particle  ejection  velocity.  Faster  parti¬ 
cles  impact  with  higher  speeds,  in  some  cases  up  to  4  km/s  (though 
still  not  as  high  as  ~20  km/s  expected  for  heliocentric  impactors). 
On  Deimos,  however,  impacts  above  1  km/s  are  relatively  rare. 
Almost  no  high-velocity  (>1  km/s)  impacts  are  noted  from  near- 
apoapsis  positions.  81%  of  impacts  are  clustered  in  the  0.4- 
0.8  km/s  region,  implying  that  low-velocity  reaccretions  to  Deimos 
are  relatively  common. 

Next,  the  relationship  between  impact  velocities  and  the  time 
to  impact  is  examined.  No  significant  acceleration  effect  with  time 
or  particular  links  to  orbital  geometry  are  noted  (Fig.  3).  However, 
the  contrast  between  the  two  bodies  is  again  evident.  On  Phobos, 
the  majority  of  impacts  occur  in  less  than  100  years;  subsequent 
impacts  become  less  frequent  as  time  increases.  This  implies  that 
the  500-year  timescale  selected  is  adequate  to  capture  the  majority 
of  Deimos-to-Phobos  material  transfer.  On  Deimos,  however, 
impacts  continue  to  build,  with  the  flux  of  impacts  remaining  rel¬ 
atively  constant  even  at  the  end  of  the  500-year  timescale.  There¬ 
fore,  it  seems  that  while  the  majority  of  mass  transfer  to  Phobos 
occurs  early  on,  reaccretions  to  Deimos  likely  continue  out  to  the 
104year  timescale  hypothesized  by  Soter  (1971).  This  also  makes 
it  likely  that  a  large  number  of  the  particles  still  flying  at  the  end 
of  the  500-year  simulation  will  end  up  reaccreting  to  Deimos.  As 
a  consequence,  it  is  unlikely  that  increasing  the  computation  time 
will  qualitatively  change  our  conclusions. 

The  flight  path  angle  (FPA)  is  the  angle  between  the  incoming 
velocity  vector  and  the  position  vector  defined  by  the  surface  of 
the  planet,  and  can  be  calculated  as  (Curtis,  2013): 


cotT  = 


esm& 

1  +  e  cos  0 


(10) 


where  e  is  the  eccentricity  of  the  impact  trajectory  and  is  the  true 
anomaly. 

The  relationship  between  impact  velocity  and  FPA  reinforces 
the  rarity  of  high-speed  reaccretion  events  on  Deimos  (Fig.  4).  On 
Phobos,  the  frequency  distribution  of  impactor  velocities  and  flight 
path  angles  suggests  that  impacts  created  by  Deimos  ejecta  can 
vary  from  oblique,  classically  secondary  impacts  to  direct  cratering 
events. 

Finally  we  test  the  fidelity  of  our  results  by  investigating  the 
likelihood  of  continuing  collisions  beyond  the  chosen  500-year 


timeframe.  Fig.  5  shows  time  curves  for  particles  still  flying  and 
particles  impacting  Phobos  or  Deimos.  As  expected,  Phobos 
impacts  taper  off  with  time,  and  the  total  number  of  Phobos 
impacts  (y)  fits  well  ( R 2  >  0.995)  to  a  logarithmic  distribution 
defined  by  y  =  178.4  In (t/t0),  where  t  is  time  in  years  and  the  time 
constant  t0  is  ~1 0.46  years.  Deimos  impacts  continue  to  increase 
and  fit  well  (R2  >  0.995)  to  a  distribution  defined  by  y  =  469 
In (t/toX  where  the  time  constant  t0  is  —47.9  years.  Interestingly, 
the  ratio  between  the  time  constants  for  Phobos  and  Deimos  is 
similar  to  the  ratio  of  their  orbital  periods. 

Reaccretions  to  Deimos  are  therefore  expected  to  continue,  but 
for  how  long?  The  shape  of  the  graph  for  particles  still  in  flight  is  in 
a  logarithmic  decrease;  when  extrapolated  (Fig.  5,  right)  it  takes 
approximately  10,000  years  for  the  number  of  particles  still  in 
flight  to  drop  below  10%  of  the  total  number  of  particles  generated. 
This  result  agrees  well  with  predictions  made  by  Soter  (1971). 
However,  a  word  of  caution  is  appropriate  here.  Soter  and  this 
study  both  neglect  effects  of  solar  radiation  pressure.  While  a 
second-order  effect  on  the  500-year  timescale,  it  can  play  a  signif¬ 
icant  role  across  longer  time  periods,  and  could  decrease  the  time 
to  impact  (Klacka,  2002).  Based  on  this  we  do  not  expect  our  con¬ 
clusions  to  change  qualitatively  with  an  increase  in  propagation 
time  and  corresponding  decrease  in  number  of  particles  still  flying. 

By  analyzing  the  Martian  system  within  the  framework  of  an 
analytical  restricted  three-body  problem,  previous  work  finds  that 
essentially  all  ejecta  from  either  Phobos  or  Deimos  will  be  reac¬ 
creted  to  the  Moon  of  origin  (Dobrovolskis  and  Burns,  1980; 
Soter,  1971 ).  Our  differing  results  suggest  that  the  four-body  ephe- 
merides  formulation  is  critical  to  full  understanding  of  the  orbital 
dynamics. 


4.  Mass  transfer  between  Martian  satellites 

To  investigate  the  geologic  impact  of  mass  transfer  between 
Phobos  and  Deimos,  we  need  to  convolve  the  probability  distribu¬ 
tions  from  Fig.  1  with  an  appropriate  mass-velocity  distribution. 
Advanced  scaling  laws  developed  from  numerical  methods  exist 
(e.g.  Leinhardt  and  Stewart,  2012),  but  given  the  uncertainties 
associated  with  several  key  parameters  we  prefer  a  more  transpar¬ 
ent  and  simpler  approach.  For  gravity-dominated  cratering  the  vol¬ 
ume  ejected  faster  than  a  given  velocity  is  (Holsapple,  1993; 
Housen  and  Holsapple,  2011): 
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Fig.  2.  (Left)  Impact  velocity  at  Phobos  and  (Right)  impact  velocity  at  Deimos  versus  particle  ejection  velocity  from  Voltaire.  While  Phobos  exhibits  several  high-velocity 
impacts,  high-velocity  impacts  at  Deimos  are  relatively  rare  and  are  primarily  clustered  below  1  km/s  (compare  to  ~20  km/s  heliocentric  impactor  velocity).  The 
discontinuity  at  700  m/s  in  both  graphs  is  due  to  an  increase  in  impacts  to  Mars  at  that  velocity  range  for  certain  Deimos-Phobos  orbital  geometries  (Fig.  1). 
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Fig.  3.  (Left)  Impact  velocity  at  Phobos  and  (Right)  impact  velocity  at  Deimos  versus  duration  of  particle  flight.  While  impacts  to  Phobos  are  frequent  in  the  first  100  years 
post  impact,  they  begin  to  taper  off  toward  the  end  of  the  examined  duration;  impacts  to  Deimos,  on  the  other  hand,  continue  at  a  relatively  constant  pace,  implying  that 
while  500  years  captures  the  bulk  of  Phobos  mass  transfer,  reaccretions  to  Deimos  will  likely  continue  to  the  104year  timescale  (Davis  et  al.,  1981;  Soter,  1971). 


where  vej  is  the  ejection  velocity  (from  Fig.  1)  and  R  is  the  final 
radius  of  Voltaire. 

A  sand-like  surface  is  well  represented  by  v  =  1.2  (Melosh, 
1989).  Experimental  results  that  determine  mass-velocity  distribu¬ 
tions  for  impacts  into  granular  targets  find  Cej  =  0.25  (Hermalyn 
and  Schultz,  2013).  Their  results  correlate  well  to  the  literature; 
Andrews  (1975),  Cintala  et  al.  (1999)  and  Stoffler  et  al.  (1975)  find 
Cej  values  between  0.25  and  0.36.  We  adopt  Cej  =  0.3  and  v  =  1.2. 
For  Phobos’  density  we  use  p  =  1.9g/cc  (Avanesov  et  al.,  1989; 


Rosenblatt  et  al.,  2008;  Schmedemann  et  al.,  2014).  From  (11), 
the  mass  ejected  within  each  velocity  bin  is  calculated  (Fig.  6).  In 
total,  3xl09kg  is  ejected  from  Deimos  between  400  and 
2000  m/s. 

Using  (11)  and  the  acceleration  due  to  gravity  of  Deimos 
(0.003  ms2),  the  total  mass  excavated  faster  than  escape  velocity 
is  6.1  x  1011  kg;  0.5%  of  this  total  is  therefore  ejected  in  the  400- 
2000  m/s  velocity  range.  Of  this  0.5%,  what  percentage  reaches 
Phobos?  From  integrating  the  mass  delivered  per  velocity  bin 
(Fig.  6)  and  averaging  it  across  impacts  at  all  GCs  (Fig.  1),  approx¬ 
imately  3.5  x  108  kg  impacts  Phobos.  This  is  12%  of  the  total  mass 
released  in  the  400-2000  m/s  range.  This  is  also  21%  of  the  mass 
not  still  in  orbit  at  the  end  of  the  simulation  (“still  flying”,  Fig.  1). 
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Fig.  4.  (Left)  Impact  FPA  at  Phobos  and  (Right)  impact  FPA  at  Deimos  versus  the  impact  velocities  on  the  respective  bodies.  Apart  from  reinforcing  the  fact  that  Deimos 
impacts  are  primarily  low-velocity,  we  also  see  a  wide  distribution  in  flight  path  angles.  There  are  several  low-velocity,  low-FPA  impacts  that  should  create  oblique  or 
secondary  crater  morphology,  and  several  high-velocity,  high-FPA  impacts  that  will  exhibit  direct  or  primary  crater  morphology. 


If  all  this  mass  were  ultimately  to  impact  Phobos,  the  total  impact¬ 
ing  mas  would  be  6  x  108  kg,  the  likely  maximum  value.  Therefore, 
12-21%  of  the  mass  that  can  reach  Phobos  does  end  up  on  Phobos 
on  a  104-year  timescale. 

The  same  analysis  yields  7.2  x  108  kg  impacting  Deimos,  which 
is  25%  of  the  total  mass  released  in  the  400-2000  m/s  range  and 
42%  of  the  mass  in  this  range  not  still  flying.  This  yields  an  upper 
bound  of  1.2  x  109  kg  impacting  Deimos.  Ejecta  launched  at  lower 
velocities  than  400  m/s  cannot  reach  Phobos  and  will  mostly  re¬ 
impact  Deimos.  Mars  only  receives  5-9%  of  ejecta  in  the  400- 
2000  m/s  range.  The  remaining  ejecta  escapes  the  Mars  system 
into  heliocentric  space. 

4.1.  Comparison  with  background  flux 

We  seek  to  place  the  sesquinary  mass  flux  into  perspective  by 
comparing  it  to  the  estimated  background  mass  flux  from  mete- 
oroidal  impacts.  The  exact  flux  of  small  meteoritic  bodies  at  Mars 
orbit  is  not  known,  so  we  derive  a  formulation  dependent  on  a 
characteristic  timescale  that  eliminates  the  Mars  mass  flux  quan¬ 
tity.  Brown  et  al.  (2002)  use  data  from  geostationary  satellites 
around  the  Earth  to  estimate  a  power  law  relationship  between 
the  number  of  objects  colliding  with  the  Earth  per  year  (N)  with 
diameters  of  at  least  D,  of  the  form: 

logN  =  c0  -  d0  logD  (12) 

where  c0  =  1.568  ±  0.03,  d0  =  2.70  ±  0.08.  Assuming  the  same  power 
law  distribution  for  Martian  system  bodies,  for  N  =  Npianet  this  can  be 
reformulated  as: 


kg  per  year.  Differentiating  (13),  substituting  into  (14)  and  integrat¬ 
ing  to  Dmax  =  D,  the  mass  accumulated  by  a  generic  planetary  body 
in  kg/yr  from  asteroidal  flux  is: 


M planet  (d  >  D)  —  Cpianet 


d0  1 

m 

[3  -  d0J 

V  6  J 

(15) 


where  d0  <  3.  Zahnle  et  al.  (2003)  derive  a  relationship  for  the 
impact  rate  of  a  satellite  compared  to  its  planet.  Applied  to  the  Mars 
system,  this  is: 

Npho(d  >  D)  =  Nmarsf  ph0  (16) 


where  D  is  the  diameter  of  the  largest  impactor  incident  to  Phobos 
and: 


f Pho 


Ri 


R-marsO-Pho 


(17) 


where  a  is  the  distance  from  Mars.  The  resulting  values  are  5  x  10-7 
and  4  x  10-6  for  Deimos  and  Phobos  respectively.  Combining  (13) 
and  (16)  for  Phobos: 


r'  _  f Pho  AfSt 

Lph0  ~  fAd^^mars 
Ust 


(18) 


where  Nsfars  signifies  Nmars(d  >  Dst).  The  assumption  here  is  that 
Stickney  is  the  largest  impact  to  have  occurred  on  Phobos  over  its 
history  and  the  diameter  of  the  Stickney-forming  impactor  is  Dst. 
Substituting  (18)  into  (15),  and  applying  (13)  to  Nsfars : 


Mpho(d  >  Dst) 


Cmarsf pj 


D 


3-d0 

st 


(19) 


Npianet (d  >  D)  =  CpianetD~d°  (13) 

where  Cpianet  =  10c°.  As  we  will  show,  the  value  of  this  constant  does 
not  matter  for  our  analysis.  Assuming  spherical  impactors  with 
diameter  D  and  density  pf,  the  incremental  number  of  impactors 
dN  per  year  results  in  a  mass  flux  increment  of 

dM  =  lpinD3dN  (14) 


Next,  we  define  a  characteristic  timescale  t  defined  such  that  r 
years  elapse  between  Voltaire-size  collisions  on  Deimos.  From  (16) 
and  the  definition  of  N: 

T  =  1/Noei  =  WDeK°L)  (20) 

where  N“jrs  is  Nmars(d  >  D„0(mire)  for  Deimos.  Expanding  this  accord¬ 
ing  to  (13): 
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Fig.  5.  Number  of  (Left)  particles  still  flying  and  (middle)  impacts  to  Phobos  and  Deimos  with  time.  The  curves  are  well  behaved,  with  no  unexpected  jumps.  Impacts  to 
Phobos  can  be  seen  to  be  tapering  off,  while  reaccretions  to  Deimos  continue  to  rise.  These  are  expected  to  continue  until  no  more  particles  are  still  flying.  (Right) 
Extrapolation  of  the  logarithmic  decrease  in  the  left  plot.  The  decrease  fits  well  ( R 2  >  0.995)  to  a  logarithmic  distribution  defined  by  y  =  -563  In (t/t0),  where  t0  ~  25,000  years. 
Using  this  curve,  it  takes  10,000  years  for  the  particles  still  flying  to  drop  below  500,  agreeing  with  Soter  (1971). 


^  —  V( f DeiCmarsDvoi  ) 


(21) 


Across  t  years,  combining  (19)  and  (21),  the  poorly  known  (and,  on 
long  timescales,  time-variable)  annual  mass  flux  delivered  to  Mars 
represented  by  Cmars  cancels.  The  total  mass  accreted  by  Phobos 
due  to  background  impacts  on  a  Voltaire  timescale  is: 


Wlpho,acc  —  T Mpho  —  r 

J  Dei 


do  1 

m 

[3  -  d0J 

V  6  J 

■'i3-d0 


(22) 


Eq.  (22)  is  only  dependent  on  the  diameter  of  the  impactors  and  the 
slope  of  the  size-frequency  distribution.  For  Phobos,  the  largest 
crater  is  Stickney,  with  a  170-m  likely  impactor  size  (Asphaug 
and  Melosh,  1993).  For  Deimos,  the  Voltaire  impactor  diameter  is 


estimated  from  gravity-dominated  scaling  relations,  rearranged 
from  Cintala  and  Grieve  (1998)  and  Schmidt  and  Housen  (1987) 
and  similar  to  Zahnle  et  al.  (2003): 


Dyol  — 


^0.862Dt 


1.2821 


(23) 


where  the  subscript  i  denotes  the  impactor  that  created  Voltaire, 
subscript  D  denotes  Deimos,  units  are  CGS  and  Dt  is  the  diameter 
of  the  Voltaire  transient  crater,  taken  to  be  1.95  km.  Asphaug  and 
Melosh  (1993)  assume  an  impact  of  3  km/s  for  the  Stickney  impact; 
assuming  the  same  impact  velocity  and  impactor  density,  (23) 
yields  a  Voltaire  impactor  diameter  of  25  m.  From  (22),  for  an 
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Fig.  6.  (Left)  Total  mass  faster  than  ejection  velocities  (Eq.  (11)).  (Right)  Mass  ejected  per  velocity  bin  in  tens  of  millions  of  kg.  Velocity  bins  correspond  to  Fig.  1. 


asteroid-type  impact  (p,  =  2.6  g/cc;  Barnhart  and  Nimmo,  201 1 ),  the 
mass  accreted  by  Phobos  from  Solar  System  impactors  between 
Voltaire-scale  impacts  is  approximately  2.9  x  109kg;  it  would  be 
less  if  the  impactor  were  assumed  to  have  the  lower  characteristic 
density  of  Phobos.  Comparing  this  to  the  3.5  x  108  kg  that  impacts 
Phobos  during  every  Voltaire  impact,  the  fraction  of  Deimos  mate¬ 
rial  delivered  to  Phobos  (F)  represents,  on  average,  0.12  of  the  mate¬ 
rial  accreted  to  Phobos  from  direct  Solar  System  impactor  flux. 

Using  a  similar  derivation  for  Deimos,  we  find  that  only 
2xl08kg  is  accreted  to  Deimos  by  Solar  System  impactors 
between  Voltaire-size  events.  Deimos  receives  less  material  than 
Phobos  because  the  focusing  factor  (Eq.(17))  is  smaller.  This  flux 
is  exceeded  greatly  by  the  reaccreted  mass  ejected  during  a 
Voltaire-size  impact;  for  ejecta  with  velocities  400-2000  m/s, 
F  =  3.6.  However  this  velocity  range  is  only  0.5%  of  the  total  mass 
thrown  into  orbit.  Ejecta  with  velocity  <  400  m/s  has  insufficient 
velocity  to  reach  Phobos  (or  Mars)  and  will  likely  reaccrete  to  Dei¬ 
mos.  Thus,  the  true  value  of  F  across  all  ejecta  is  likely  ~700  for 
Deimos.  In  other  words,  on  Deimos,  mass  reaccreted  during  a 
Voltaire-size  impact  greatly  exceeds  mass  naturally  accreted  from 
Solar  System  impactors.  However,  sesquinaries  from  a  large  impact 
on  Deimos  provide  only  a  minor  contribution  to  the  flux  at  Phobos. 

As  a  reality  check,  we  calculate  the  applicability  of  (12)  to  Mars 
system  impacts.  The  frequency  of  a  Voltaire-sized  impact  is  z  =  1  / 
Ndei  years,  or  134  Ma.  We  discuss  the  issue  of  the  apparent  age  of 
Voltaire  further  below.  Using  (21),  the  ratio  of  the  timescale  for  a 
Voltaire-forming  impact  on  Deimos  to  a  Stickney-forming  impact 
on  Phobos  is  0.049.  Therefore,  the  frequency  of  a  Stickney-sized 
impact  is  approximately  every  2.7  Ga,  which  seems  quite 
reasonable. 

As  an  additional  check,  we  estimate  the  mass  flux  to  Mars  from 
studies  in  the  literature.  Chappaz  et  al.  (2011)  find  that  the  mass 
flux  from  Mars  to  Phobos  can  be  estimated  at  0.25  pg/m3/yr  or 
0.2  kg/yr  across  the  inner  moon’s  surface  area  assuming  a  mixing 
depth  of  0.5  m.  The  mass  flux  of  Solar  System  projectiles  to  Phobos 
is  numerically  found  to  be  k  =  40-2400  times  greater  than  the  flux 
from  Mars  ejecta,  with  a  Monte  Carlo  preferred  value  of  k  =  195 
(Ramsley  and  Head,  2013).  The  preferred  value  yields  a  mass  flux 
at  Phobos  of  38.5  kg/yr  from  asteroids,  comets  and  meteoroids. 
At  this  rate,  across  134  Ma,  Phobos  accumulates  5.2  x  109  kg  from 


direct  impacts,  within  a  factor  of  two  of  the  2.9  x  109  kg  derived 
above. 

While  independent  of  Mars’  meteoroidal  flux,  the  results  are 
admittedly  susceptible  to  the  assumed  size  of  the  Voltaire  impac¬ 
tor.  The  size  of  Stickney  implies  that  the  gravity  regime  approxima¬ 
tion  is  likely  appropriate.  However  the  short  timescale  derived  for 
the  Voltaire  impact  implies  that  it  may  be  deeper  in  the  strength 
regime  than  initially  assumed.  If  true,  the  Voltaire  impactor  would 
have  to  be  larger  than  25  m,  which  would  reduce  quantitative  esti¬ 
mates  of  F,  the  fraction  of  Deimos’  mass  flux  with  respect  to  the 
Solar  System  impactor  mass  flux. 

Assuming  that  the  impactor  was  incident  at  a  45°  angle 
(Holsapple,  1993;  Melosh,  1989),  we  can  estimate  the  diameter 
of  an  impactor  for  strength-dominated  cratering  by  rearranging 
Eq.  (12)  from  Zahnle  et  al.  (2008)  to  yield: 

Dvo,  =  1.027  (24) 

where  Y  represents  the  dynamic  strength  of  the  body;  other  vari¬ 
ables  are  as  in  (23).  Melosh  (1989)  uses  a  value  of  Y=2  MPa,  the 
observed  yield  stress  at  crater  collapse,  to  estimate  the  gravity / 
strength  transition  on  the  Earth  and  the  Moon.  Adopting  this  value 
in  (24),  the  diameter  of  the  Voltaire  impactor  is  88  m.  However,  the 
value  of  Y  is  uncertain;  for  the  low-density  Deimos,  it  is  unlikely 
that  the  yield  strength  is  as  high  as  2  MPa.  For  an  order  of  magni¬ 
tude  change  in  Y  ranging  from  0.2  to  20  MPa,  the  diameter  of  the 
impactor  varies  from  40  to  190  m. 

Is  strength  or  gravity  scaling  more  appropriate  for  modeling  the 
Voltaire  impact?  Gravity  can  be  a  factor  on  Solar  System  bodies  as 
small  as  400  m  (Love  and  Ahrens,  1996).  Phobos  and  Deimos  have 
similar  compositions,  bulk  densities  and  accelerations  due  to  grav¬ 
ity  (Davis  et  al.,  1981;  Szeto,  1983),  and  several  arguments  for 
modeling  cratering  on  Phobos  in  the  gravity-dominated  regime 
are  detailed  in  Asphaug  and  Melosh  (1993)  and  Asphaug  et  al. 
(2015a).  Finally,  while  some  authors  have  used  the  wide  distribu¬ 
tion  of  ejecta  on  Deimos  to  surmise  that  strength-scaling  may  be 
appropriate  for  Deimos  (Lee  et  al.,  1986),  we  suggest  that  this  glo¬ 
bal  distribution  could  instead  be  a  function  of  the  large  amount  of 
mass  reaccreted  over  500  yr  timescales  (Section  3).  However,  even 
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if  the  Voltaire  impact  is  in  the  strength  regime,  this  would  not 
change  our  qualitative  results,  i.e.,  (1)  that  the  sesquinary  mass 
transfer  is  a  relatively  small  fraction  of  meteoroidal  impacts  to 
Phobos  and  (2)  that  the  mass  reaccreted  to  Deimos  from  large  Dei- 
mos  impacts  exceeds  the  meteoroidal  mass  flux  to  Deimos.  For 
instance,  increasing  the  diameter  of  the  Voltaire  impactor  to 
88  m,  the  qualitative  results  become  F=  0.004  for  Phobos  and 
F  =  24  for  Deimos  across  all  ejecta  velocities.  Therefore,  uncertainty 
on  where  the  Voltaire  impact  lies  in  the  gravity/strength  regime 
does  not  affect  our  qualitative  conclusions. 

5.  Discussion 

5.1.  Importance  of  inter-moon  mass  transfer  and  re-accretion 

The  central  result  of  this  work  concerns  the  relative  mass  trans¬ 
fer  during  impacts  on  Deimos.  We  have  found  that  a  Voltaire-sized 
impact  on  Deimos  does  transfer  mass  from  Deimos  to  Phobos,  with 
sufficient  velocity  to  create  primary  crater  morphology,  discussed 
further  below.  When  viewed  with  reference  to  the  Solar  System 
impactor  flux,  the  sesquinary  mass  transfer  is  not  significant,  and 
is  likely  in  the  10%  range.  However,  compared  to  the  25  ppm  of 
Mars  material  in  Phobos  regolith  estimated  by  Chappaz  et  al. 
(2011),  one  out  of  ten  particles  originating  from  Deimos  is  a  large 
number  and  is  of  interest  to  Phobos  lander  mission  concepts  (e.g. 
Udrea  et  al.,  2015,  2016). 

While  the  placement  of  Voltaire  within  the  strength-gravity 
domain  can  cause  some  uncertainty  in  impactor  size,  we  have 
shown  that  the  total  ejecta  mass  reaccreted  to  Deimos  is  likely  to 
exceed  the  background  mass  flux  to  Deimos  (F>20)  greatly.  A 
Voltaire-sized  impact  could  therefore  feasibly  resurface  large  parts 
of  the  Moon,  erasing  the  previous  geological  record.  Dating  the  sur¬ 
face  of  Deimos  may  be  more  challenging  than  previously  sus¬ 
pected;  the  surface  age  may  better  represent  the  age  of  either 
Voltaire  or  the  similarly  sized  Swift  crater.  Further,  an  11 -km  con¬ 
cavity  on  the  southern  end  of  Deimos  is  hypothesized  to  be  a  pos¬ 
sible  impact  scar  from  an  ancient,  very  large  collision  (Lee  et  al., 
1986;  Thomas,  1993;  Thomas  et  al.,  1996).  If  true,  this  would  have 
resulted  in  the  transfer  of  significant  sesquinary  mass  transfer  to 
Phobos,  and  a  complete  resurfacing  of  Deimos’s  surface. 

5.2.  The  spectral  dichotomy  of  Phobos 

Phobos  exhibits  two  distinct  spectral  units,  one  of  “redder”  ori¬ 
gin  and  one  of  “bluer”  origin,  a  distinction  that  likely  stems  from  a 
compositional  difference  (Lee,  2009;  Rivkin  et  al.,  2002).  It  has 
been  proposed  that  the  red  unit  is  a  wide-spread  shallow  layer 
superimposed  on  a  blue  base  that  is  perhaps  more  representative 
of  Phobos  composition  at  depth  (Murchie  and  Erard,  1996; 
Murchie  et  al.,  2008).  Possible  mechanisms  for  the  superimposition 
of  red  material  are  briefly  outlined  in  Section  1  (Britt  and  Pieters, 
1988).  These  include  the  hypothesis  that  the  spectrally  red 
“veneer”  may  be  ejecta  accreted  from  Deimos  rather  than  Mars, 
as  suggested  by,  e.g.,  Smith  et  al.  (2015). 

Our  results  for  the  distribution  of  low-velocity,  oblique-angle 
impacts  (Fig.  4)  support  the  existence  of  trajectories  that  could, 
in  theory,  deposit  a  “veneer”  of  red  Deimos  material  across  Pho- 
bos’s  surface.  However,  the  inter-moon  mass  flux  is  relatively  small 
compared  to  the  Solar  System  impactor  flux,  which  likely  has  a 
greater  effect  on  the  global  surface  geology,  particularly  in  the 
100+  Ma  since  the  last  Voltaire-sized  impact.  Therefore  we  believe 
it  unlikely  that  the  red  veneer  of  Phobos  is  of  Deimos  origin.  Recent 
spectral  analysis  by  Thomas  et  al.  (201 1 )  supports  this  finding  with 
evidence  of  subsequent  impacts  penetrating  the  blue  unit  near 
Stickney  to  reveal  redder  material.  This  observation  casts  doubt 


on  the  idea  that  the  blue  unit  may  be  representative  of  Phobos  at 
depth  (Basilevsky  et  al.,  2014).  Ultimately,  sample  return  from  both 
Phobos  and  Deimos  will  conclusively  establish  which  spectral  unit 
is  representative  of  depth;  our  results  suggest  that  the  surface  unit 
is  unlikely  to  originate  from  Deimos. 

5.3.  Primary  versus  secondary  impact  morphology 

Due  to  the  dearth  of  classic  secondary  impact  features  such  as 
radial  crater  chains  or  herringbones,  it  has  been  concluded  that 
few,  if  any  Phobos  craters  are  secondary  in  origin  (Thomas  and 
Veverka,  1980b;  Thomas,  1979).  A  similar  conclusion  was  reached 
for  Deimos.  A  limit  of  10  m/s  on  maximum  re-impact  velocity  was 
proposed  by  Thomas  (1998)  and  continues  to  be  used  in  analysis  of 
Phobos’  geology  (Murray  et  al.,  2006;  Schmedemann  et  al.,  2014; 
Smith  et  al.,  2015).  However,  every  sesquinary  impact  studied  here 
occurs  at  speeds  above  100  m/s.  Though  still  orders  of  magnitude 
below  heliocentric  impactor  velocities,  Fig.  4  shows  several  high- 
velocity,  high-FPA  particles  that  could  create  craters  indistinguish¬ 
able  from  primary  impact  craters  on  Deimos,  and  particularly  on 
Phobos. 

On  Deimos,  over  80%  of  impacts  cluster  in  the  0.6  ±  0.2  km/s 
region,  implying  that  the  majority  of  reaccretions  are  low  velocity 
(subsonic).  However  for  Phobos  sesquinary  particles  can  arrive 
with  either  subsonic  or  supersonic  impact  velocities.  One  would 
therefore  predict  a  wide  range  in  the  resulting  crater  morphology; 
a  comprehensive  image  survey  of  Phobos  may  be  able  to  distin¬ 
guish  between  the  different  crater  morphologies.  We  find  several 
low-velocity,  low-FPA  particles  incident  to  Phobos  that  could  cre¬ 
ate  classic  oblique  or  chained  secondary  impact  morphology;  a 
Phobos  image  survey  by  Smith  et  al.  (2015)  finds  several  craters 
and  deposits  likely  originating  from  such  low-velocity  impacts. 

The  escape  speed  from  Mars  at  Phobos’  orbit  is  3  km/s,  which  is 
exceeded  by  the  fastest  sesquinary  impacts  (Fig.  2).  These  impacts 
can  create  ejecta  of  their  own  that  may  subsequently  be  lost  from 
the  Mars  system.  Mass  loss  from  Phobos  is  one  possible  explana¬ 
tion  for  why  outlines  of  ejecta  blankets  are  not  conspicuous  on 
Phobos  (Lee  et  al.,  1986).  Additional  simulations  of  Phobos- 
centered  ejecta  dynamics  are  needed  to  confirm  this,  and  are 
planned  as  future  work. 

5.4.  Deimos  surface  smoothness ,  brightness  and  sesquinary  impact 
gardening 

Very  little  ejecta  escapes  large  bodies  (e.g.  Earth),  with  the 
majority  redeposited  locally  as  a  continuous  ejecta  blanket.  Most 
ejecta  escape  from  very  small  bodies  (e.g.  Phoebe  Burns  et  al., 
1996),  never  to  be  reaccreted.  Ejecta  dynamics  on  Deimos  present 
an  interesting  bridge  between  these  two  regimes,  with  sesquinary 
effects  appearing  to  be  important.  Ejecta  escapes  but  is  then  reac¬ 
creted  on  a  timescale  of  up  to  hundreds  of  years,  resulting  in  a  glo¬ 
bal,  near-isotropic  redistribution  of  sesquinary  ejecta. 

For  a  Voltaire-sized  impact,  almost  all  of  the  ejecta  material 
launched  at  velocities  <400  m/s  will  ultimately  reaccrete  to 
Deimos  (Soter,  1972;  Thomas  et  al.,  1996).  From  Eq.  (11)  above  this 
represents  6.1  x  1011  kg,  or  about  0.5  m  thickness  of  material 
distributed  evenly  over  the  entire  surface.  We  estimate  that  the 
Voltaire  and  Swift  impacts  together  could  have  added  on  the  order 
of  1  m  of  fresh  regolith  or  crater  fill  to  Deimos.  It  has  been 
suggested  that  the  smoother  and  brighter  surface  of  Deimos  is 
due  to  crater  fill  of  5-7  m  depth  (Thomas,  1979,  1993;  Thomas 
et  al.,  1996;  Veverka,  1978).  Given  our  estimate,  we  suggest  that 
reaccreting  sesquinary  mass  provides  at  least  a  partial  explanation 
for  the  origin  of  this  crater  fill  material. 

The  impact  velocity  distributions  (Fig.  2)  also  suggest  a  tertiary 
ejecta  effect.  Though  most  sesquinary  impacts  to  Deimos  are 
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relatively  low  speed  (<1  km/s),  these  are  still  significantly  higher 
than  the  escape  velocity  of  Deimos  (5.5  m/s,  Table  1)  and  could 
potentially  launch  additional  ejecta  in  their  turn  (see  below).  Sim¬ 
ilarly,  ejecta  from  Deimos  can  impact  Phobos  with  enough  mass 
and  speed  to  create  craters  and  excavate  Phobos  mass,  which  could 
then  enter  Mars  orbit. 

On  Deimos,  sesquinary  impacts  represent  a  large  mass  flux  rel¬ 
ative  to  the  background  flux  and  have  ~km/s  impact  velocities 
(Fig.  2).  The  result  is  likely  to  be  the  production  of  further  subor¬ 
bital,  ballistically  emplaced  ejecta.  This  mechanism  could  con¬ 
tribute  to  the  smooth  appearance  of  the  Deimos  surface.  Energy 
and  momentum  transfer  from  reaccretion  impacts  might  even  set 
off  downslope  movement  noted  on  Deimos  (Thomas  and 
Veverka,  1980a),  though  admittedly  the  efficiency  of  the  impacts 
at  initiating  this  process  is  unknown  and  should  be  constrained 
in  the  future. 

The  possibility  of  sesquinary  impact  gardening  makes  the  evo¬ 
lution  of  regolith  in  between  major  collision  events  a  complex  pro¬ 
cess  on  both  Martian  moons.  Determining  the  exact  nature  of 
Deimos  material  mixed  with  Phobos  regolith  is  one  of  the  primary 
science  objectives  of  lander  concepts  in  development  for  Phobos 
(Udrea  et  al.,  2016,  2015).  The  methods  applied  in  this  work  are 
further  applicable  to  understanding  regolith  development  and  dust 
belts  on  small  bodies  within  planetary  gravitational  wells,  such  as 
Phoebe  and  Iapetus  in  the  Saturnian  system. 
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extends  through  the  intersection  between  the  planet’s  equator 
and  its  prime  meridian,  with  the  y-axis  completing  the  right  hand 
system. 

Fig.  7  illustrates  the  relationship  between  the  PCI  and  PCPF 
frames.  PCI  can  be  rotated  into  the  PCPF  frame  around  the  z-axis 
with  the  rotation  matrix: 


PCPFftPCI 


cosw  sinw  0 
-sinw  cosw  0 
0  0  1 


A3.  Topocentric  Horizon  frame  (also  SEZ:  South-East-Zenith  frame) 

The  topocentric  horizon  frame  is  adapted  from  the  South- 
East-Zenith  (SEZ)  frame  as  defined  by  Vallado  (2013);  we  refer  to 
the  two  interchangeably  here.  The  definition  of  the  Topocentric 
Horizontal  frame  assumes  a  sphere  centered  at  the  center  of  mass 
of  Deimos  and  tangent  to  the  origin  of  the  frame,  which  is  the  cen¬ 
ter  of  the  impact  site  (Voltaire).  The  x-axis  is  aligned  with  the 
meridian  that  passes  through  the  center  of  Voltaire  and  points 
south.  The  y-axis  is  defined  such  that  the  x-y  plane  is  tangent  to 
the  surface  of  Deimos  at  the  center  of  Voltaire  and  points  along 
the  local  latitude  circle.  Completing  the  right-handed  system,  the 
z-axis  points  radially  outward  from  the  impact  site  towards  the 
“local”  zenith.  The  local  horizon  forms  the  fundamental  plane  for 
this  system,  i.e.,  the  plane  defined  by  the  south  and  east  axes.  It 
should  be  noted  that  there  is  a  subtle  difference  between  the 
definition  of  the  impact  site’s  latitude  by  geodetic  or  astronomical 
standards  (Vallado,  2013);  these  become  identical  by  imposing  the 
assumption  of  a  perfectly  spherical  impacted  body  (Deimos).  The 
low  gravitational  acceleration  at  the  surface  of  Deimos  has  a 
negligible  effect  on  the  speed  of  the  ejecta,  so  the  assumption  of 
a  uniform  spherical  geometry  is  justified. 

Fig.  8  illustrates  the  relationship  between  SEZ  and  PCPF.  sej  is 
the  elevation  angle  of  the  ejection  velocity  vector  from  the  hori¬ 
zontal,  defined  as  0  <  sej  ^  90°.  From  the  Z-model  formulation 
(Eq.  (1 )),  sej  =  35.4°.  pej  is  the  azimuth  of  the  ejection  velocity  vector 
and  is  measured  from  the  North,  positive  clockwise  as  viewed  from 
above  the  site  such  that  0  s:  Pej  360°,  and  is  sampled  at  11  posi¬ 
tions  across  this  range  separated  by  30°.  The  streamlines  defined  in 
Z-model  frame  can  be  rotated  into  SEZ  using  the  relationships: 

Vsouth  =  Vr  COS (180°  fiej) 


Appendix  A.  Coordinate  transformations 

This  section  details  coordinate  transformations  necessary  to  use 
surface-centered  Z-model  streamlines  in  the  Mars-centered  simu¬ 
lation.  Three  planet-based  frames  are  used  and  detailed  here. 

A  2.  Planet-Centered  Inertial  (PCI)  frame 

The  origin  of  the  PCI  frame  is  the  center  of  the  planet:  Mars, 
Deimos  or  Phobos.  The  positive  x-axis  points  toward  the  vernal 
equinox,  the  positive  z-axis  extends  through  the  North  Pole  of 
the  planet  and  the  y-axis  completes  the  right  hand  system.  In  this 
definition,  the  North  Pole  is  that  pole  of  rotation  that  lies  on  the 
north  side  of  the  invariable  plane  of  the  Solar  System  (Archinal 
et  al.,  2010).  The  planetary  system  model  described  in  Section  2.3 
is  placed  with  reference  to  the  Mars  Centered  Inertial  (MCI)  frame. 

A2.  Planet-Centered  Planet-Fixed  (PCPF)  frame 

Like  the  PCI  frame,  the  origin  of  the  PCPF  frame  is  also  the  center 
of  the  planet,  and  shares  its  z-axis  definition.  However,  the  x-axis 


VEast  =  Vr  sin(180°  -  Pej) 


Fig.  7.  Illustration  of  the  relationship  between  the  Planet-Centered  Inertial  (PCI) 
and  Planet-Centered  Planet-Fixed  (PCPF)  frames. 
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Fig.  8.  Illustration  of  the  relationship  between  the  Planet-Centered  Planet  Fixed 
(PCPF)  frame  and  the  South-East-Zenith  (SEZ)  frames.  The  fundamental  plane  is 
highlighted  (blue).  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 

where  vr  is  defined  by  Eq.  (2).  Subsequently  two  rotation  matrices, 
the  first  around  the  y-axis  and  the  second  around  the  z-axis,  are 
required  to  rotate  the  SEZ  frame  into  the  PCPF  frame  (specifically, 
the  Voltaire  SEZ  frame  into  the  Deimos-Centered  Deimos-Fixed 
frame).  The  rotation  matrices  are: 


"  cos  Ai 

-  sin  Ai 

O' 

sin^j 

0 

cos  cpt ~ 

DCDFffSEZ  _ 

sin  Ai 

cos  Ai 

0 

0 

1 

0 

0 

0 

1 

-  cos  cpt 

0 

sin  cpt 

Finally,  we  can  convert  the  DCDF  velocity  coordinates  to  MCI. 
By  manipulating  the  basic  kinematic  equation  for  the  position  vec¬ 
tor  of  an  ejected  particle  in  the  DCDF  frame,  rPCDF,  in  terms  of  the 
known  z)PCDF  we  get  the  expression: 

vf'  =  +  iP™  +  Ma(0DCDF  X  r?CDf 

^  rjfa  =  v%CI,  where  v%CI  is  the  velocity  vector  of  Deimos  in  the  MCI 
frame.  This  is  determined  from  Deimos  ephemerides.  MCIc oDCDF  is  the 
angular  rate  of  the  DCDF  frame  with  respect  to  MCI;  this  is  the  rota¬ 
tion  rate  of  Deimos.  rfCDF  is  the  position  of  the  impact  site  in  DCDF 
coordinates. 

Following  a  similar  process,  for  the  impact  of  ejecta  with  Pho- 
bos  the  process  is  reversed  to  obtain  the  velocity  vector  at  impact 
in  the  Phobos  Centered  Phobos-Fixed  (PhCPhF)  frame.  With  similar 
notation  as  used  above  this  relation  is: 

^PhCPhF  _  pMCI  _  fiMCI  _  MCI  ^PhCPhF  x  jiPhCPhF 

The  validity  of  these  equations  has  been  checked  with  MICE,  a  com¬ 
mercial  level  interface  created  by  JPL/Caltech  to  SPICE  ephemeris 
information  from  NASA’s  Navigation  and  Ancillary  Information 
Facility  (naif.jpl.nasa.gov)  (Acton  et  al.,  2002). 
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